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We discuss the phenomenon of the smooth dynamical gravity induced crossing of the phantom 
divide line in a framework of simple cosmological models where it appears to occur rather naturally, 
provided the potential of the unique scalar field has some kind of cusp. The behavior of cosmological 
trajectories in the vicinity of the cusp is studied in some detail and a simple mechanical analogy is 
presented. The phenomenon of certain complementarity between the smoothness of the spacetime 
geometry and matter equations of motion is elucidated. We introduce a network of cosmological 
histories and qualitatively describe some of its properties. 

PACS numbers: 98.80.Cq, 98.80.Jk 



I. INTRODUCTION 

The discovery of cosmic acceleration pj has stimulated 
a construction of a class of dark energy models 0- [J [| 
describing this effect. This dark energy should possess a 
negative pressure such that the relation between pressure 
and energy density w is less than —1/3. Some observa- 
tions indicate that the present day value of the param- 
eter w < — 1 provides the best fit. The corresponding 
dark energy has been named phantom dark energy 
According to some authors, the analysis of observations, 
permits to specify the existence of the moment when the 
universe changes the value of theparameter w from that 
the region w > —1 to w < —1 y|, 13 ■ This transition is 
called "the crossing of the phantom divide line" . 

It is easy to see, that the standard minimally coupled 
scalar field cannot give rise to the phantom dark energy, 
because in this model the absolute value of energy den- 
sity is always greater than that of pressure, i.e. \w\ < 1. 
A possible way out of this situation is the consideration 
of the scalar field models with the negative kinetic term. 
Thus, the important problem arising in connection with 
the phantom energy is the crossing of the phantom di- 
vide line. The general belief is that while this crossing 
is not admissible in simple minimally coupled models its 
explanation requires more complicated models such as 
multifield ones or models with non-minimal coupling be- 
tween scalar field and gravity (see e.g. H3). 

In our preceding paper [TlT ] we have described the phe- 
nomenon of the change of sign of the kinetic term of the 
scalar field implied by Einstein equations. Now we try 
to answer the question how general is this phenomenon. 
In other words, is it some curiosity arising due to a very 
particular choice of the form of a potential of a scalar 
field and of initial conditions or it is rather a typical phe- 
nomenon. In doing this, we introduce a notion of the net- 
work of cosmological histories and qualitatively describe 
its properties. Admittedly, our research looks purely aca- 



demical, but we stress that the persisting signals in favor 
of the dark energy equation of state parameter w < — 1 
;9j justifies the interest in this topic and, in particular, 
in our systematic approach. The structure of the paper 
is as follows: Sec. II is devoted to a brief recapitulation 
of the technique of reconstruction of scalar (phantom) 
potentials and the transformation of the unique scalar- 
phantom field (i.e. the crossing of the phantom divide 
line; in Sec. Ill we discuss the role of initial conditions 
for the phantom divide line crossing phenomenon; Sec. 
IV is devoted to the analysis of a simple mechanical sys- 
tem with the potential with a cusp, which is strictly mo- 
tivated by the requirement of the possibility of crossing 
of the phantom divide line within a framework which al- 
lows a unique scalar field; in Sec. V we introduce the 
notion of network of cosmological histories; the section 
VI contains concluding remarks In the Appendix A we 
shall present simple solvable examples of such networks. 
In the Appendix B we shall describe some potentials as 
functions of time or scalar field, for which exact solution 
of the Einstein equations is possible. 



II. RECONSTRUCTION OF SCALAR 
(PHANTOM) POTENTIALS AND DYNAMICAL 
CROSSING OF THE PHANTOM DIVIDE LINE 

Let us remember the main points of our approach. We 
have considered the minimally coupled scalar field with 
an exponential potential 



V {4>) ~ exp 



(1) 



Such potentials are widely studied in cosmology 0] and 
have a particular solution corresponding to a special 
choice of initial conditions, which implies the power-law 
behavior of the cosmological radius 



a(i) = aot^ 1 ^, 



(2) 



2 



or, in other terms, 



h(t) 



3(l + w)t' 



(3) 



where h — d/a is the Hubble variable. Changing the ini- 
tial conditions (i.e. the value of </> corresponding to some 
value of <fi) one shall have other laws of cosmological evo- 
lution, which cannot be presented explicitly, but which 
have the same qualitative behavior as @. 

However, one can consider the potential |^ not as a 
function of the scalar field 4>, but as a function of the 
cosmic time parameter t as was done in [TlL [l3L [T^ . fT5| . In 
this case the Einstein (Fricdmann) equation is equivalent 
to the linear second-order differential equation 



for the volume function 

i;(t) = a 3 (t). 



(4) 



(5) 



The potential V(t) corresponding to the exponential po- 
tential JIJ and to the choice of the initial conditions pro- 
viding the evolution (J5J) is 



Vi(t) 



2(1 - w) 



(6) 



9(1 + w) 2 t 2 ' 

Equation Q with the potential JHJ has a general solution 



ip(t) = ^ii Ql 
2 



■■02*° 



a 2 



1 + iu' 

w — 1 



(7) 



where ipi and tp2 are nonnegative constants. If '02 = 0, 
one reproduces the cosmological evolution while if 
ipi = the cosmological evolution requires the pres- 
ence of the phantom scalar field with the negative sign 
of the kinetic term. The most interesting situation re- 
veals when, both the coefficients ipi and ^2 are nonzero. 
Simple calculation shows that in this case the Hubble 
variable has a maximum point where, the kinetic term is 
forced to change its sign. 1 . 

Now, fixing the relation between these coefficients one 
can ask what is the corresponding form of the potential 
as a function of the scalar field 4>. The explicit form of 
this function cannot be found, but we have studied [TH 
its asymptotic behavior at singularities and in the point 



to, where the kinetic term changes the sign and (de)- 
phantomization occurs, or, in other words, where the 
universe crosses the phantom divide line. In the neigh- 
borhood of the moment to, which depend on the relation 
between the coefficients ip\ and ip2 UM > the potential has 
the form 



v{4>) 



2(1 



H = 



9(1 + w) 2 [t + ( 
V8(l-w)(3-™) 2 0i02io T 




(1 



(8) 



The scalar field behaves correspondingly as 
0~(i-M 3/2 - 



(9) 



III. INITIAL CONDITIONS AND PHANTOM 
DIVIDE LINE CROSSING 

In the previous section we have implicitly chosen such 
initial conditions so that simultaneously the scalar field 
and its time derivatives vanish at the point to, which 
in turn implies the crossing of the phantom divide line. 
However, let us fix the form of the potential as a function 
of the scalar field 4> and see how the behavior of the 
universe will change depending on the initial conditions. 
Suppose that when the cosmic parameter is close to to 
the potential behaves as 



V{4>) 



(A + Bcj) 2 / 3 ) 2 ■ 
The scalar field can have the following form: 

<p(t) = 4>o + Mt--to) a , 



(10) 



(11) 



where a > 1 to provide the vanishing of the time deriva- 
tive of the scalar field at the point to- The situation 
described above and corresponding to the (de) - phan- 
tomization of the scalar field 01 realizes when <po = 
and a = 3/2. Let us consider the cases when (f>o =/= 0, 
analyzing the Klein-Gordon equation 



+ 3h<t> - 



dV 



= 0. 



The terms here behave as 

= -to)"" 2 , 



(12) 



(13) 



1 Notice that in papers 111! the technique of the reconstruction of 
the potentials as functions of the redshift parameter was devel- 
oped. The reconstruction approach is (to a large extent) model 
independent. So, it may show that the phantom boundary is 
crossed in nature, but do not specify a class of potentials as 
functions of <j> in which this can be described in a consistent way 



4> = a4>i(t - t ) 



dV 4B 



a-l 



(14) 



(15) 



3 



h=Jz(x 2 <t>l(t-to) 2a - 2 



(A- 



(16) 



Let us consider the case a > 2. In this case the first 
and the second terms of Eq. (|12[1 vanish at t — > while 
the term (|15|) is a constant. Thus, we have come to a 
contradiction and this case should be disregarded. 

Now, considering the case a = 2, we see that while 
the term containing the first derivative of the scalar field 
vanishes the term with the second derivative becomes a 
constant (2<f>i) to be equated to the the term (|15(l up to 
the sign. Therefore, we obtain the relation 



2B 



3^1/3(^4 + 5^2/3)3 



(17) 



It is easy to see that in this case the time derivative of 
the Hubble parameter h behaves as 



{t-tf 



(18) 



in contrast to the behavior h ~ t — to, corresponding to 
the case 4>q = 0, considered in Apparently, the be- 
havior (|18fl means that the Hubble variable has at the 
moment to & n inflection point instead of the extremum 
point and hence, the kinetic term for the scalar field con- 
serves its sign. 

In the case when 1 < a < 2 the second derivative term 
diverges while other terms are finite, so this behavior 
cannot be realized. 

Until now we have considered the cosmological trajec- 
tories having at some moment the vanishing time deriva- 
tive <fi of the scalar field, while <f> could have zero or 
nonzero values. However, there is also another class of 
trajectories characterized by the fact that at some mo- 
ment t\ the scalar field vanishes, while its first time 
derivative is different from zero. Such trajectories are 
characterized by the equation: 



0(t) =V 1 (t-t 1 )+V 5/3 (t-t l ) 



5/3 



(19) 



where the constants v\ and D5/3 satisfy the condition 



V5/3 



6B 



5A 3 ' 



1/3 • 



(20) 



which can be easily obtained from Eqs. Q12[l and (|15fl . 
Calculating the second time derivative of the Hubble vari- 
able, one can see that it behaves as h(t) ~ (t — ti) -1 / 3 
that means that higher curvature invariants including the 
curvature derivatives diverge and one encounters the so 
called soft singularity. The problem of removal of such 
singularities is of interest by intself, but will not be con- 
sidered here. 

Thus, let us summarize the results of the above consid- 
eration. We have found three types of the trajectories. 
For the first family of the trajectories the time derivative 



of the scalar field vanishes at the moment to ■ This family 
is a one-parameter one and is parametrized by the value 
of the scalar field (f>o , which should be different from zero 
at the moment to- All the nonzero values of (f>o require 
a = 2, the value of 4>\ is fixed by relation (|17f) . the Hub- 
ble variable has an inflection point and the kinetic term 
conserves its sign. 

The second family contains the only trajectory with 
4>o = 0. Now the term (|15(l becomes singular and be- 
haves as (t — to)~ a / 3 , the term with the first derivative 
vanishes, while the second derivative term behaves like 
(t — to) a ~ 2 ■ Consistency requires a = 3/2 and the value 
of 4>\ is uniquely fixed. The Hubble variable has at the 
moment to an extremum and crossing of the phantom 
divide line occurs. The third family is characterized by 
Eqs. lfT§|) . (|2U|> and suffers from a soft cosmological sin- 
gularity. 

It seems that in some sense one can consider the pos- 
sibility of the crossing of the phantom divide line as was 
described in 0] as exceptional and the corresponding 
initial conditions as having measure zero. Nevertheless, 
we would like to show that considering broader set of 
potentials and initial conditions, one can come to the 
conclusion, that conditions providing the crossing are, 
roughly speaking, commensurable to those excluding it. 
We postpone the corresponding consideration until the 
fifth section, while the next section will be devoted to an 
analysis of the simple mechanical model of the particle 
moving in the potential with a cusp. 



IV. MECHANICAL ANALOG: A PARTICLE 
MOVING IN A POTENTIAL WITH A CUSP 

Let us consider a one-dimensional problem of a classi- 
cal point particle moving in the potential 



V(x) = 



(1+ £2/3)2' 

where Vb > 0. The equation of motion is 

4U 



3(1 + X 2 /3)3 a; l/3 



0. 



(21) 



(22) 



We consider three classes of possible motions character- 
ized by the value of the energy E. The first class consists 
of the motions when E < Vo- Apparently, the particle 

with x < 0, x > or with x > 0, x < cannot reach the 

3/2 

point x — and stops at the points =F I y ^ — 1^ 

respectively. This class of motions corresponds to the 
class of cosmological evolutions described in the preced- 
ing chapter when the Hubble variable has an inflection 
point and the universe does not cross the phantom divide 
line. 

The second class includes the trajectories when E > 
Vo- In this case the particle crosses the point x — 
with nonvanishing velocity and this case correspond to 



4 



the cosmological evolution given by Eq. I|19|l . but natu- 
rally there is nothing singular in it, because in our simple 
mechanical problem, we do not have nothing similar to 
the Hubble variable and the spacetime geometry defined 
by its behavior. 

If we have a fine tuning such that E = Vo, we encounter 
an exceptional case. Now the trajectory satisfying Eq. 
(|22|l in the vicinity of the point x = can behave as 



x = C(t - t) 



3/2 



where 



C = ± 



16Vp 
9 



3/4 



(23) 



(24) 



and t < to. It is easy to see that independently of the sign 
of C in Eq. H24|) the signs of the particle coordinate x 
and of its velocity x are opposite and hence, the particle 
can arrive in finite time to the point of the cusp of the 
potential x = 0. 

Another solution reads as 



= C{t-t a f\ 



(25) 



where t > to- This solution describes the particle going 
away from the point x = 0. Thus, we can combine the 
branches of the solutions ill'MI) and (|25() in four different 
manners and there is no way to choose if the particle ar- 
riving to the point x = should go back or should pass 
the cusp of the potential (|21|l . It can stop at the top 
as well. Such a "degenerate" behavior of the particle in 
this third case is connected with the fact that this trajec- 
tory is the separatrix between two one-parameter fami- 
lies described above. At the moment there is not yet any 
strict analogy between this separatrix and the cosmologi- 
cal evolution describing the phantom divide line. In order 
to establish a closer analogy and to understand what is 
the crucial difference between mechanical consideration 
and general relativistic one, we can try to introduce a 
friction term into the Newton equation (|22|l 



4% 



jx - 



3(1 + 3^/3)3^1/3 



= 0. 



(26) 



It is easy to check that if the friction coefficient 7 is a 
constant one does not have a qualitative change in re- 
spect to the discussion above. Let us asuume for 7 the 
dependence 



then 



and 



7 = — x 
' 2 



7 = — 3ii. 



(27) 



(28) 



(29) 



The trajectory arriving to the cusp with vanishing veloc- 
ity is still described by the solution (|23|l . Consider the 
particle coming to the cusp from the left (C < 0. It is 
easy to see that the value of gamma at the moment to 
tends to zero, while its second derivative 7 given by Eq. 
(EH is 



> 0. 



(30) 



Thus, it looks like the friction coefficient 7 reaches its 
minimum value at t = to. Let us suppose now that the 
particle is coming back to the left from the cusp and its 
motion is described by Eq. i|25|) with negative C. A 
simple check shows that in this case 



7(*o) 



-C 2 < 0. 



(31) 



Thus, from the point of view of the subsequent evolution 
this point looks as a maximum for the function 7(i). In 
fact, it means simply that the second derivative of the 
friction coefficient has a jump at the point t = to- It is 
easy to check that if instead of choosing the motion to the 
left, we shall move forward our particle to the right from 
the cusp (C > 0), the sign of gamma(t ) remains nega- 
tive as in Eq. (|31|) and hence we have the jump of this 
second derivative again. If one would like to avoid this 
jump, one should try to change the sign in Eq. I|29|l . To 
implement it in a self-consistent way one can substitute 
Eq. (123 by 



7 = 3 



~ + V{x) 



and Eq. HJJl by 



Wo 



- 72; - 



3(1 + x 2/3)3 a ,l/; 



= 0. 



(32) 



(33) 



In fact, it is exactly that what happens automatically 
in cosmology, when we change the sign of the kinetic en- 
ergy term for the scalar field, crossing the phantom divide 
line. Naturally, in cosmology the role of 7 is played by 
the Hubble variable h. The jump of the second derivative 
of the friction coefficient 7 corresponds to the divergence 
of the third time derivative of the Hubble variable, which 
represents some kind of soft cosmological singularity (this 
singularity is even softer than that considered in the pre- 
ceding section for the family of cosmologies characterized 
by Eq. l|19fl where already the second time derivative of 
the Hubble variable was divergent). 

Thus, one seems to confront the problem of choosing 
between two alternatives: 1) to encounter a weak singu- 
larity in the spacetime geometry; 2) to change the sign 
of the kinetic term for matter field. In this paper we 
pursue the second alternative insofar as we privilege the 
smoothness of spacetime geometry and consider equa- 
tions of motion for matter as less fundamental than the 
Einstein equations (see also more detailed discussion in 



5 



V. POTENTIAL AS A COSMIC TIME 
FUNCTION AND THE NETWORK OF 
COSMOLOGICAL HISTORIES 



Indeed, all the previous considerations were based on 
the treatment of the potential as a fixed function of the 
scalar field and changing initial conditions for this scalar 
field. As was already emphasized 0,0, E| one can use 
an alternative treatment of the potential V as a function 
of time. In the example considered in |Tlj , it was shown 
that for V(t) instead the initial conditions which exclude 
the crossing of the phantom divide line are exceptional. 

Let us try to imagine a set of possible cosmologies as 
a two-dimensional surface (see the picture below) where 
the vertical lines represent the potentials as functions of <j> 
while the horizontal lines represent the potentials as func- 




Vo{<l>) 



Network of cosmological histories 

tions of t. 

We shall call it network of cosmological histories. Mov- 
ing along these lines corresponds to changing of initial 
conditions keeping fixed the form of the potential. For 
example, for the vertical line corresponding to the expo- 
nential potential JIJ (see the line Vo{4>) in the picture, the 
black circle A corresponds to the power-law cosmologi- 
cal evolution given by formulae J5J) and QJ) there are no 
initial conditions allowing the crossing of the phantom di- 
vide line. Now, fixing the initial conditions providing the 
cosmological evolution @, @ we obtain the potential 
V\(t) ~ 1/t 2 and moving along the corresponding hori- 
zontal line one finds that at all the points, excluding the 
point of intersection A, there is the (de)-phantomization 
effect. Let us stop now at some point of the horizon- 
tal line (say, B) and fix the corresponding form of the 
potential V^(^) as a function of cj>. Then, one can begin 
moving along the vertical line, corresponding to this form 
of V2 (</>). The intersection point B corresponds to such 
choice of initial conditions when <fr and phi vanish simul- 
taneously at some moment to. And this is the only point 
of the vertical line, where the crossing of the phantom 
divide line occurs as was shown above. To come to this 
conclusion we have used only some general properties of 
the form of the potential around a candidate moment of 



the crossing of the phantom divide line. 

Now, we undertake a similar analysis showing that for 
a horizontal line (potential as a function of t) almost 
all the initial conditions imply the crossing effect. Let us 
consider a potential V3 (t) corresponding to the horizontal 
line intersecting the vertical line V2 (</>) at the point C. We 
shall look for a cosmological evolution where the Hubble 
variable has an extremum at the moment to. A regular 
potential can be represented around this moment as 



V(t) = ^(U + Ui(t-to)). 



Then Eq. acquires the form 



= (U + U lT )fl>, 



(34) 



(35) 



where r = t — to- We shall look for two independent 
solutions in the vicinity of r = 0; in the form 



(36) 



n=0 



The recurrent relation between the coefficients a n is 
(n + 2)(n + l)a n+2 = U a n + Uia n -\. 



(37) 



We define the first independent solution by fixing of the 
first two coefficients as do = 1, a\ = 1 while for the sec- 
ond solution we choose ao = \,a\ — (one can easily 
check that these solutions are independent, calculating 
the first terms of their Wronskian) . The leading terms of 
the general solution of Eq. (|35|l has the following form: 



i>(T) = ipi 1 + r + 



U t 2 ([/ + £/i)r 3 



G 



+1p2 1 + 



U T 2 



UlT 3 

6 



(38) 



where ipi and 1P2 are nonnegative. We shall need also the 
expressions for the first and second time derivatives of 
the solution l|38|) : 



tp = tpi 1 + V t + 



{Vo+Vxy 



+ V>2 V T + 



(39) 



t/j = ijj! (U + (U + U x )r) + ^2 (Uo + U x t). (40) 
The time derivative of the Hubble variable has the form 



• 1 ■i/jijj — tp" 2 
3 i> 2 



(41) 



It is enough to study the numerator of expression (1411) 

faj>-# 2 = {1>i + ik) 2 Uo-ift 



(42) 



6 



To have at the moment r = the extremum of the 
function h(r) its time derivative should behave as h ~ r, 
that means that Eq. (|42|l implies that 

W>i+^2) 2 E/o = Vi\ (43) 

and 

t/i ^ 0. (44) 

The condition Q43|l determines the relation between 
the coefficients ipi and tp2 in terms of the coefficient Uq. 
Thus, for a potential V(t) fixing the moment of time to 
we fix the constants Uq and U\. If the condition (I44|) 
is satisfied, choosing the constants tpi and ip2 satisfying 
condition l|43l) one obtains a cosmological evolution un- 
dergoing the crossing of the phantom divide line. 

VI. CONCLUDING REMARKS 

Thus, we have seen that considering the potential as a 
function of the scalar field having a cusp we should choose 
exceptional initial conditions to describe a smooth cos- 
mological evolution undergoing the phantom divide line 
crossing. We have seen also that there is some kind of 
complementarity between the smoothness of the space- 
time geometry and that of the structure of matter action. 

On the hand, considering the potential as a function 
of time we encounter just the following situation: initial 
conditions implying the crossing of the phantom divide 
line are typical. That means that roughly speaking in the 
set of all possible cosmological trajectories those crossing 
and non-crossing the phantom divide line have measures 
of the same order. We are not able to construct a rigorous 
description of this set of cosmological histories, but the 
treatment undertaken above permits to make some qual- 
itative remarks. First of all let us notice that a network 
of vertical lines representing the potential as a function 
of the scalar field and horizontal lines , representing po- 
tentials as functions of the cosmic time, is not closed. 
Namely, let us leave a vertical line at the point A and 
make some walk along the horizontal line to the point 
B. At the point B we shall have a cosmology with the 
crossing. Then let us make some walk along the vertical 
line from the point B to some point C. The cosmol- 
ogy corresponding to the point C is free from crossing. 
Now, traveling along the horizontal line from point C we 
shall always have trajectories with crossing and, hence, 
we shall not have an opportunity to return to the ini- 
tial vertical line, which should have cosmologies without 
crossing. (The empty circle D in our picture signifies the 
absence of intersection between the respective vertical 
and horizontal lines). 

There is another interesting question. It is well known 
that for the exponential potential all the initial con- 
ditions imply the cosmologies free from the crossing phe- 
nomenon. The explanation of this fact is a very simple 
one: the effect of smooth crossing can have place 



only for potentials which at some point have a diver- 
gent derivative with respect to the scalar field, while the 
derivative of the exponential potential is always regular. 

In conclusion we would like to say that the above inves- 
tigation shows that the smooth crossing of the phantom 
divide line in simple cosmological models could be con- 
sidered as a rather natural phenomenon more than an 
exotic anomaly. The network of cosmological histories 
constructed for the illustration of this fact seems to be 
interesting object on its own. The points of this network 
represent functions h(t) on which two transformations 
generated by means of scalar field potentials, treated as 
functions of 4> or of t. Graphically they were represented 
as displacements along vertical and horizontal lines. The 
network discussed in this paper was constructed begin- 
ning from the simple cosmological evolution @. Notice 
that our initial point of network |j2J) is characterized by 
some value of the equation of state parameter w and is 
nothing but a power-law evolution The points ly- 
ing on a horizontal line intersecting this initial power-law 
point (the point A on our picture), are described explic- 
itly [llj and have more complicated structure. As far as 
other points of the network are concerned we have stud- 
ied above some of their local properties, but we are not 
able to describe their global topology. It seems unlikely 
that two different power-law points belong to the same 
network. If this statement is confirmed one has got a 
family of networks characterized by different values of 
the initial constant equation of state parameter w. In 
the Appendix A we give an explicit illustration of the 
simplest network containing an exponential expansion, 
while in the Appendix B we give some other examples of 
potentials for which it is possible to find exact solutions. 
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Appendix A. Examples of simple networks of 
cosmological histories 

We present here a particularly simple and in some way 
degenerate example of the network of cosmic histories 
which is discrete and contains only three points. Let us 
consider de Sitter universe with 

h(t) = H . (45) 

The corresponding potential as is well-known and as can 
be easily obtained from the equation (see e.g. ^lj and 
references therein) 

V=\ + h\ (46) 
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has the form 

v(4>) = Hi 

while the initial condition is 

6 = 0. 



(47) 



(48) 



This initial condition can be defined at any moment of 
time and the evolution shall conserve it, because the 
Klein-Gordon equation l|12|l for the constant potential 
has the form 



6 + 3h6 = 0. 



(49) 



Notice that in this case the Klein-Gordon equation has 
the same form for usual scalar and phantom. Now, treat- 
ing the potential as a constant function of time we can 
easily resolve for this case Eq. Q: 



tp = ipi exp(3-ff i) + i>2 exp(-3H t). 
The Hubble function correspondingly is 



(50) 



tpi exp(3//ot) + Wi cxp(—3Hot) 
while its time derivative is 



h(t) 



(V>i exp(3H t) + V2 exp(-3#o*)) 2 



(52) 



As usual the constants tpi and are nonnegative. Cor- 
respondingly we have three opportunities. If ip2 — this 
is our initial point describing infinitely expanding flat de 
Sitter universe. If ipi = 0, one has an infinitely con- 
tracting flat de Sitter universe. If both the coefficients 
are positive, the Hubble variable is growing and this sit- 
uation could be realized only in the universe where the 
scalar field has a negative kinetic term (phantom). It is 
easy to understand that a change of the relation 4>2 / ipi is 
simply equivalent to the shift of the cosmic time variable 
t. Thus, we have only one cosmic history described by 
phantom field with the constant potential. In this history 
the universe begins its evolution at t = — 00 with a quasi- 
de Sitter contraction. It arrives to a minimum possible 
value of the cosmological radius at some moment t m i n . 
This moment can be easily found from Eq. (|51|l however, 
its value does not have a particular meaning just like the 
relation tp2 /ipi ■ After the moment t m i n begins an expan- 
sion which at t — > 00 becomes quasi-de Sitter. Making a 
convenient shift of the time variable we can write down 
Eq. I|51|l in a specially simple form: 



h(t) = H Q t&nh3H t. 



(53) 



Thus, instead of a horizontal line we obtain two discrete 
points. 

Now, let us treat the constant potential 147|) as a 
function of 6 changing the initial condition for the time 
derivative of this field. It is easy to show that for 6 7^ 



there two cosmological histories given by the following 
expression for the Hubble variable: 



h(t) = ±H Q coth3H t. 



(54) 



These solutions describe a universe which begins its evo- 
lution from cosmological singularity at t = and then 
is infinitely expanding arriving to quasi-de Sitter phase 
at t — > 00 and a universe which begins its evolution at 
t = — 00 and then contracts arriving to the Big Crunch 
type cosmological singularity at t = 0. It is important 
to notice that the choice of the initial value of the time 
derivative of the scalar field is not important, because in 
the process of evolution this value runs between and 
±00. The sign of this time derivative is not important 
either because the potential of this field is constant. 

Thus, we have presented a schematic example of an ex- 
actly solvable cosmological network, which contains only 
five discrete points. Two of these points represent an 
infinitely expanding and infinitely contracting de Sitter 
universes, two of them describe universe evolving from 
the cosmological singularity to an infinite expansion or 
vice versa, driven by a scalar field with constant poten- 
tial. The last point represents a universe which passes 
from contraction to expansion and is driven by a phan- 
tom field. All the network is characterized by an absolute 
value of the constant Hq. Indeed, beginning from — Hq 
we arrive to the same network, which was constructed 
from Hq. 

As a last remark we can give two examples of even 
more trivial networks of cosmological histories. First 
of them contains only one point. If one has a nega- 
tive constant potential V = —Hq, than the only cos- 
mological history present in the network is that given by 
h(t) = — iJotan3i?oi describing the evolution from the 
Big Bang at t = —ti/QHq to the Big Crunch at t = tt/QHq 
passing through the point of maximal expansion at t = 0. 
The second network consists of three points. Consider- 
ing massless scalar field with vanishing potential, one can 
see that there are three opportunities: if the time deriva- 
tive of the scalar field is nonzero, then one has either a 
universe which begins its evolution from the Big Bang 
singularity and expands infinitely or a universe which 
contracts finishing in the Big Crunch singularity. If the 
time derivative of the scalar field is zero, one has a static 
Minkowski universe. 



Appendix B. Some solvable potentials 

The concept of the network of cosmological histories 
introduced in this paper is based on the simultaneous 
use of the notion of potential as function of time and 
as function of the scalar field. Unfortunately, it is much 
more difficult to find exactly solvable examples for poten- 
tials as functions of scalar fields than for those, treated 
as functions of time. Thus, we, first give another exam- 
ple of the potential as function of 6 (see, for additional 
details p^l. 
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Let us consider the cosmological evolution given by the 
Hubble function 



and 



h(t) = iJ coth 



3H (1 +w)t 



while the cosmological radius behaves as 

/. u 3ff (l + W )^ 2/t3(1+ - )1 
a(t) = a sum ) 



(55) 



(56) 



This evolution occurs in the universe filled with two per- 
fect fluids: one of them is a cosmological constant and 
the other is a barotropic fluid with the equation of state 
parameter w. The scalar field potential realizing this 
evolution has the form flTl : 



V{<j>) =Hl[l 



1 — w 



■ sinh 



3v / TTw</> 



2 2 
The potential as a function of time is 

{l-w)Hl 



V(t) = H$ + 



2 sinh 



2 3H (l+w)t 



(57) 



(58) 



The general solution of Eq. J2J with potential l|58|) can 
be written as usual as 



+?/>2 sinh 



3H {l + w)t 



2/(1+™) 



/^(sinh 3g ° (1 2 +w) 



,x -4/(1+™) 



(59) 



The explicit form of the solution J7J) is very simple for 
some special values of the parameter w. We shall consider 
the case w = 0. Now, the solution is 

tl>(t) = ipi sinh 2 — I" "02 ^coth sinh3i?oi 

(60) 

If -02 = 0, we come back to our starting point: a cosmo- 
logical evolution (|53|l . (|5*oT) for w = 0. 

The case when both the coefficients ipi and ^2 are dif- 
ferent from zero is rather cumbersome, while its qual- 
itative behavior can be understood studying the case 

^i = o,v^o. 

Let us rewrite the solution l|60l for this case in the form 



-oo < t < -- 



3H 



sh2. 



The Hubble variable can be expressed as 

= 3ff (l- y + y 2 )sign(t) 
v / y 2_ T(2-y) ' 

and its time derivative reads 

9ff 2 (2/ 3 + 3y 2 -6y+l) 



h(t) 



where 



(2/2-1) {2-yf 



y = cosh3H t. 



(63) 



(64) 



(65) 



(66) 



One can see by simple algebraic treatment that the func- 
tion h(t) can have the only zero at some moment to in the 
interval (|62|l . i.e. when < y < 2. Such a point describes 
a transition from a phantom-type contraction to a non- 
phantom contraction ending in a Big Crunch. The evolu- 
tion in the interval l|63(l describes instead a non-phantom 
contraction beginning at t = — oo in a de Sitter stage and 
ending at t = —-^j^arccosh2 in the Big Crunch. 

Analogously, one can show that in the case ip2 < 0, one 
can have either non-phantom expansion beginning from 
the Big Bang and ending in an infinite de Sitter uni- 
verse or the transition from the non-phantom expansion 
to phantom expansion culminating in a Big Rip singular- 
ity at t = jjj^arccosh2. The inclusion of the term with 
V>1 7^ can change some features of the cosmological evo- 
lution. However, the qualitative and numerical analysis 
of the corresponding differential equations shows that the 
number of the moments of time when the universe crosses 
the phantom divide line does not change. 

For example, for the case when both the coefficients 
tpi and ip2 are positive and ipi = 3V>2 the evolution is 
presented in Fig. 1. Here, the universe begins its evolu- 



h[t] 



Sinh [3H t/2] 



cos h Mo* 

ip(t) = 1P2 ^(2 - cosh3ff i). (61) 

sinh — if- 

Let us consider first the case ip2 > 0. Then, the volume 
function ip(t) will be nonnegative in the intervals: 

< t < -}—arccosh2, (62) 
3-Ho 



FIG. 1: h(t) dependence at = 3 

tion at t = in the phase of phantom-type contraction, 
which at some moment (when h(t) crosses the line h = 0) 
is substituted by the phase of phantom-type expansion. 
After some period, the universe crosses the phantom di- 
vide line (when h(t) has its maximum value) and enters 
into the phase of an infinite normal expansion. 
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As we have already mentioned above, the number of 
the scalar field potentials treated as functions of the 
scalar field, for which it is possible to find an exact so- 
lution is very limited. Instead, one can find some poten- 
tials treated as functions of the cosmic time parameter, 
for which one can find complete exact solutions. It is 
connected with the linearity of the corresponding equa- 
tion (0J , which is nothing but a Schrodinger equation for 
a system with zero energy. Thus, we can borrow some 
solvable examples of such potentials from non-relativistic 
quantum mechanics (see e.g. 0,0, US)- We shall give 
here a couple of examples of such potentials. For the 
potential 



V(t) 



the solution of Eq. (@J looks like 
'A" 



ijj = ip t exp 



■^it exp 



.4 



(67) 



(68) 



For the potential 



V(t) = 



B 1 



1 



(l + < 2 ) 5 



(69) 



the solution reads 



ip — ipi\/l + t 2 exp(Barctant) 
+ip2 yl + t 2 exp(— Sarctant) . 



(70) 



The investigation of cosmologies based on the solutions 
(|68|l and l|7U|) might represent some interest but it is out- 
side of the scope of the present paper. 
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